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Abstract 

We study a class of spherically symmetric Stephani cosmological models in the presence of a 
• self-interacting scalar field in both classical and quantum domains. We discuss the construction 

of 'canonical' wave packets resulting from the solutions of a class of Wheeler- De Witt equations 
in the Stephani Universe. We suggest appropriate initial conditions which result in wave packets 
^q' containing some desirable properties, most importantly good classical and quantum correspondence. 

We also study the situation from de-Broglie Bohm interpretation of quantum mechanics to recover 
the notion of time and compare the classical and Bohmian results. We exhibit that the usage of the 
canonical prescription and appropriate choices of expansion coefficients result in the suppression of 
the quantum potential and coincidence between classical and Bohmian results. We show that, in 
some cases, contrary to Friedmann-Robertson- Walker case, the bound state solutions also exist for 
all positive values of the cosmological constant. 
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In recent years observations show that the expansion of the Universe is accelerating in the present epoch 
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PQ contrary to Friedmann-Robertson- Walker (FRW) cosmological models, with non-relativistic matter 
and radiation. Some different physical scenarios using exotic form of matter have been suggested to 
resolve this problem [2j [3l SI [6l [7] . In fact the presence of exotic matter is not necessary to drive 
an accelerated expansion. Instead we can relax the assumption of the homogeneity of space, leaving 
the isotropy with respect to one point. The most general class of non-static, perfect fluid solutions of 
Einstein's equations that are conformally flat is known as the "Stephani Universe" [HI [H] • This model 
can be embedded in a five-dimensional flat pseudo-Euclidean space [HI EH E] and its three dimensional 
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spatial sections are homogeneous and isotropic [12] . Recently, quantum spherically symmetric Stephani 
cosmological models in the presence of the perfect fluid have been studied in Refs. [121 Q3] . In these 
works the Schutz's variational formalism |15H16j is applied to recover the notion of time and investigate 
the singularity avoidance at the quantum level. 

The question of construction and interpretation of wave packets in quantum cosmology and its 
connection with classical cosmology has been attracting much attention in recent years. Moreover, 
numerous studies have been done to obtain a quantum theory for gravity and to understand its connection 
with classical physics. 

In quantum cosmology, in analogy with ordinary quantum mechanics, one is generally concerned with 
the construction of wave functions by the superposition of the 'energy eigenstates' which would peak 
around the classical trajectories in configuration space, whenever such classical-quantum correspondence 
is mandated by the nature of the problem. However, contrary to ordinary quantum mechanics, a 
parameter describing time is absent in quantum cosmology. Therefore, the initial conditions would 
have to be expressed in terms of an intrinsic time parameter, which in the case of the Wheeler-DeWitt 
(WDW) equation could be taken as the local scale factor for the three geometry [17]. Also, since the 
sign of the kinetic term for the scale factor is negative, a formulation of the Cauchy problem for the 
WDW equation is possible. The existence of such a sign is one of the exclusive features of gravity with 
many other interesting implications. 

The construction of wave packets resulting from the solutions of the WDW equation has been a 
common feature of some research works in quantum cosmology [HI QlJl [20l Ell [22]. In particular, in 
references [3T] [521 123] the construction of wave packets in a Friedmann Universe is presented in detail 
and appropriate boundary conditions are motivated. Generally speaking, one of the aims of these 
investigations has been to find wave packets whose probability distributions coincide with the classical 
paths obtained in classical cosmology. In these works, the authors usually consider theories in which 
a self interacting scalar field is coupled to gravity in a Robertson- Walker type Universe. The solutions 
are obtained such that the following desirable properties are satisfied. There should be a good classical- 
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quantum correspondence, which means that the wave packet should centered around the classical path, 
the crest of the wave packet should follow as closely as possible the classical path, and to each distinct 
classical path there should correspond a wave packet with the above properties. Recently, a general 
prescription has been suggested by Gousheh et al [22] for constructing the 'canonical' wave packets 
which contain all above desired properties. They showed that there always exists a "canonical initial 
slope" (CIS) for a given initial wave function, which optimizes some desirable properties of the resulting 
wave packet, most importantly good classical-quantum correspondence. 

In this paper we deal with the subject of "initial condition" which is an important problem in 
quantum cosmology. In fact, in classical cosmology we can uniquely determine the classical initial 
conditions subject to the zero energy condition. But in quantum cosmology since the underlying equation 
(WDW equation) is a hyperbolic differential equation, we are free to choose the initial wave function 
and the initial derivative of the wave function by choosing arbitrary expansion coefficients. These 
quantities (distributions) correspond to classical initial position and initial momentum, respectively. 
Therefore, although WDW equation allows us to use different choices of initial conditions upon choosing 
different expansion coefficients, these wave functions correspond to different classical situations. This 
also happens whenever WDW-like equation appears in other theories such as varying speed of light 
quantum cosmological models |23j . Hence, a legitimate question which arises is how we can construct 
a specific wave packet which completely corresponds to its unique classical counterpart? One possible 
solution is removing the arbitrariness of the expansion coefficients and defining a certain relation between 
them. In our previous investigations [221 123] we discovered that given a particular choice of initial wave 
function, certain coefficients remain undetermined, and if we set the functional form of those coefficient 
to be the same as the determined ones, we obtain excellent classical and quantum correspondence. 

Here, we are interested to study the Stephani Universe in the presence of a scalar field. First, we 
write the reduced action near symmetry center (r ss 0) and find the corresponding hamiltonian. Then 
we obtain the Einstein's equations and WDW equations in minisuperspace. These equations can be 
solved numerically with appropriate initial conditions. In particular, we use Spectral Method (SM) [24] 
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as an accurate and stable numerical method for solving the quantum cosmology case which can be cast 
in the form of a hyperbolic PDE. The form of the scalar potential is chosen to contain some desirable 
properties like the cosmological constant, positive mass term in Taylor expansion, and to be bounded 
from below. Then, we construct the wave packets for various functional forms of the spatial curvature 
through the canonical prescription [52] . 

The paper is organized as follows: In SecJH we outline the main problem which is a case of Stephani 
cosmology where the matter is taken to be a particular type of self-interacting scalar field. We derive 
the main equations both for the classical cosmology and the quantum cosmology. We begin Sec. [3] with 
a description of the Spectral Method [24] which is a robust numerical method, and then we review the 
general prescription of canonical wave packets. We then consider various cases and solve them in both 
classical and quantum cosmological domains. In Sec. [4] we find the corresponding Bohmian trajectories 
and compare the classical and quantum solutions. In Sec.[5J we draw some final conclusions. 

2 The model 

Let us start from the Einstein-Hilbert action plus a scalar field as 

S=\l d 4 x^jU + 2 / d 3 xVhh ab K ab + f <£xj—g ( -i(V0) 2 - U{4>) J , (1) 

2 JM JdM JM V 2 / 

where K ab is the extrinsic curvature and h ab is the induced metric over the three-dimensional spatial 
hypersurface, which is the boundary dM of the four dimensional manifold M in units where 8nG = 1 
[55] , The last term of ([T]) represents the scalar field contribution to the total action. 

The metric in spherically symmetric Stephani Universe [8] [9] [12l [TO] [26] [27] has the following form 

ds 2 = ~N 2 (r, t)dt 2 + -|^T \dr 2 + r 2 {d9 2 + sin 2 6 # 2 )1 , (2) 
V z (r,t) 

where 
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is the lapse function and the functions V(r,t) and F(t) are defined as 



V(r,t) 
F(t) 



l + -k(t)r 2 , 
R(t) 



(4) 
(5) 



where k, R, and C are arbitrary functions of time [28] [29] ■ Here, k(t) plays the role of the spatial 
curvature and R(t) is the Stephani version of the FRW scale factor. Although k(t) is an arbitrary 
function of time in the Stephani model, assuming a power law relation between R(t) and k{t) makes the 
model solvable and is in agreement with the accelerating expansion of the Universe [T5irTlll2"TIl2"9"ll3T)II31| . 
However, some authors have used some thermodynamics relations to obtain a power law relation between 
these two variables |31j . Though in the spherical symmetric inhomogeneous models, the scalar field <f> 
depends on both r and t, we can consider it as an only function of time near the symmetry center r ~ 
which means = <j>{t) [3"2"l |3"3"1 151] . 

By substituting the Stephani metric J2J) in the action JT]) and choosing the curvature function k(t) 
in the form [T3J QH 123 EH EDI [31] 

k(t)=0W(t), (6) 
after dropping the surface terms and with due attention to the form of the lapse function ([3]), the final 
reduced action near r ~ takes the form 



S = J dt 

Now choosing the gauge N = 1 



k 2 R 



N 



+ 3(3NR 1+1 + NR 3 Q^ 2 - U{4>) 



we have the following Lagrangian 
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L = -3R 2 R + 3/3i? 1+7 + R 3 ( -0 2 - U(<j>) 



(7) 



(8) 



Therefore, in this limit, the Stephani Universe is equivalent to the FRW model where the curvature term 
can be chosen as an arbitrary function of time. 

The Einstein's equations for r ~ resulting from above Lagrangian with the zero energy condition 
can be written as 



R 
R 



(9) 



2 (l) + (l) +^ + ^ 2 = + (io) 

D CiTJ 

* + 3-* + ^ . 0, (11) 

where dot represents differentiation with respect to time. We require the potential U (0) to have natural 
characteristics for small 0, so that we may identify the coefficient of ^</> 2 in its Taylor expansion as a 
positive mass squared m 2 , and U(Q) as a cosmological constant A. An interesting choice of U(<p) with 
three free parameters is [l"9l [21] |2"2[ |35 [ |36 [ [37] 

U(4>) =A+^r sinh 2 (a0) + sinh(2a<£). (12) 

In the above expression m 2 = d 2 U/d(f> 2 \ ( p=o is a mass squared parameter and b is a coupling constant. 
We need to choose a 2 = g in order to separate the variables in the Lagrangian. This potential is bounded 
from below and as we shall see, prevents us from usual problem of factor ordering. Moreover, since this 
type of potential also has been used for FRW cosmological models, we can compare our solutions with 
the previous FRW results [2T1I22]. 

The Lagrangian ([5]) can be cast into a simple form using the transformations X = R 3 / 2 cosh(a0) 
and Y = i? 3 / 2 sinh(a0), which transform the term i? 3 U((j)) into a quadratic form. Upon using a second 
transformation to eliminate the coupling term in the quadratic form, we arrive at new variables u and 
v, which arc linear combinations of X and Y 



coshf? sinh6> \ / X 
sinhfl cosh6> I [ Y 



(13) 



where 



e = -tanh- 1 (^) . (14) 
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L{u, v) 



I -2 2 2\ 1-2 2 2\ ^ / 2 2\(t + 1)/3 

(U — U X U ) — (V — U) 2 V ) — -^P(u — V ) w " 



(15) 



where u>\ 2 = — 3A/4 + m 2 /2 =F V m 4 — 46 2 /2. The resulting Einstein's equations are 



u + u\u + M( 7 + i) u ( u a _ w 2)(7-2)/3 = 0i ( i 6) 



y + w f„ + ^(7 + l)v{u 2 - v 2 )^' 3 = 0, (17) 

it 2 + lo\u 2 -v 2 - cu 2 2 v 2 + 9 -l3(u 2 - v 2 )^' 3 = 0. (18) 

Equations (fl6|) and (fTT)) are the dynamical equations and (fT8|) is the zero energy constraint. The 
non-linearity of these equations for 7 7^ 2 is now apparent. The corresponding quantum cosmology is 
described by the Wheeler-DeWitt equation written as 

W{u, v) = + |1 + ulu 2 u 2 2 v 2 + \p(u 2 v 2 )^ 3 } tf( tt , v) = 0, (19) 

which arises from the zero energy condition psp . In general, this equation is not exactly solvable and 
we should resort to a numerical method [24 . 

3 Solutions for the quantum cosmology cases 

We start this section by a discussion of the numerical method that we shall use and then we outline the 
general prescription for finding CIS. The general hyperbolic PDE that we want to solve is 

where f'(u,v) is an arbitrary function. It is notable that such equations may represent a wave-like 
equation whose solution may rapidly oscillate. In such cases, the usual spatial integration routines such 
as Finite Difference Methods fail to produce a reasonable solution. Therefore, it is of prime importance 
to use a reliable, efficient and accurate numerical method [23] . 

SM [35] consists of first choosing a complete orthonormal set of eigenstates of a preferably relevant 
hermitian operator to construct the solution. Since the whole set of the complete basis has usually 
infinite elements, we make the approximation of representing the solution by only a finite superposition 
of the basis functions. By substituting this approximate solution into the differential equation, a matrix 
equation is obtained. The expansion coefficients of these approximate solutions could be determined by 
the eigenfunctions of this matrix. In this method, the accuracy of the solution is increased by choosing 
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a larger set of basis functions. Having resorted to a numerical method, it is worth setting up a more 
general problem defined by the following WDW equation 



Hip(u, v) 



d 2 d 2 



(20) 



As mentioned before, any complete orthonormal set can be used. Here we use the Fourier series basis by 
restricting the configuration space to a finite square region of sides 2L. This means that we can expand 
the solution as 



z 



mmu \ ( nirv 



L J 9j \ L J' 



i,j=l m,n 



where 



( mjru \ _ / 2 • / miru \ 

{ mnu \ 
\ L ) 



92 m 



and R r< 



COS 



1, m ^ 

2, m = 



(21) 



(22) 



By referring to the WDW equation (j2"0|) , we realize that in the Fourier basis it is appropriate to introduce 
/' as 



f'(u, v) = f(u, v) + uj\u 2 - io\v 2 



(23) 



We can make the following expansion 



f'(u, V)tp(u, V) = ^2^2 B 'm, n ,i,j 9* {~J~) 9] (-p) 

i,j m,n 



(24) 



where B' m n i j are coefficients that can be determined once f'(u, v) is specified. By substituting (|21j24j) 
in ([20]) and using the independence of gi ( 22 p t )s and gj ( 2 p i )s we obtain 



~Z~) 



nir\' 

T) 



Am,n,i,j + B mni j — 0, 



(25) 



where 



B' 



E 



L i-L 



-L J-L 



'mnu\ /nirv\ / tti'ttu 

9* [ —r~ I 9 j ~r f (u, v)g v —— ) g jt 



L ) 



L 



L 



dudv 



Am' ,n' ,z' ,j f 



/ j KJ m,n,i,2,m' ,n' ,i' ,j' -"m' ,n',t' ,j' ■ 



(26) 
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Therefore we can rewrite (|2"5)l as 

m,n,i,j + ^ f C mni j m t n , j/ ji A m i j7l r y ji = 0. (27) 

Now, we select 47V 2 basis functions, that is m and n run from 1 to TV. It is obvious that the presence 
of the operator f'(u,v) leads to nonzero coefficients C' m n i j m , nl i , ,j, in ([27]) , which in principle could 
couple all of the matrix elements of A. Then we replace the square matrix A with a column vector A' 
with (27V) 2 elements, so that any element of A corresponds to one element of A' . This transforms (f2"7| 
to 

DA' = 0. (28) 

Matrix D is a square matrix with (27V) 2 x (27V) 2 elements which can be obtained from (j2"T|) . Equation 
(|28|) can be looked as an eigenvalue equation, i.e. DA' a = aA' a with (27V) 2 eigenvectors. However, for 
constructing the acceptable wave functions, i.e. the ones satisfying the WDW equation (f2"0)) , we only 
require eigenvectors which span the null space of the matrix D. That is, due to (|27[) we will have exactly 
27V null eigenvectors which will be linear combination of our original eigenfunctions introduced in (|2ip. 
After finding the 27V eigenvectors of D with zero eig envalue, i.e. A' k (k = 1, 2, 3, 27V), we can find 
the corresponding elements of matrix A, A^ ni ■. Therefore, the wave function can be expanded as 

where A fe s are arbitrary complex coefficients which can be fixed by the initial conditions. 

We are free to adjust two parameters: 27V, the number of basis elements and 2L, the length of the 
spatial region. This length should be preferably larger than spatial spreading of all the sought after wave 
functions. However, if 2L is chosen to be too large we loose overall accuracy. Therefore, it is important 
to note that for each TV, L should be properly adjusted [35]. 

Now to determine X k s we need to apply the initial conditions. As a mathematical point of view, 
since the underling differential equation is second order, X k s are arbitrary and independent coefficients. 
On the other hand, if we are interested in constructing the wave packets which simulate the classical 
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behavior with known classical positions and velocities, these coefficients will not be all independent 
yet. To address this issue, let us study the problem near the solution's boundary (v — 0). We can 
approximate (fT9|) near the v — 0, so up to the first order in v we have 

(>1 + ^ + w y + V* + i)/ 3 L ( ^ ) = 0i (30) 



du 2 dv 2 1 4 
This PDE is separable in u and v variables, so we can write 



ip(u,v) = ip J (u)x(v). (31) 



By substituting tp(u,v) in (|3H|) , two ODEs can be derived 

d 2 Xn{v) 



dv 2 



+ E nXn (v) = 0, (32) 



d 2 vl(u) , / , , 2 9 



/3« 2(7+1)/3 )^(«) = E nV l{u), 



ihjl +^(u 2 + -Pu^>^)^(u) = E n <pl(u), (33) 

where E n s are separation constants. These equations are Schrddinger-like equations with E n s as their 
'energy' levels. Equation ([32]) is exactly solvable with plane wave solution as 

X n (v) = a n cos {^/e^ vj + ip„ sin ( \f~E~ n v) , (34) 



where a n and (3 n are arbitrary complex numbers. Equation (|33| does not seem to be exactly solvable and 
we resort to a numerical technique. As mentioned before, SM can be used to find the bound state energy 
levels (E n ) and the corresponding wave functions (tp n (u)) with high accuracy. The general solution to 
the (|3"D|) can be written as 

ip(u,v)= 2J {A n cas{\/E^v)+iB n sm(y/E^v)^'l{u)+ {C n cos{^/E^ l v)+iD n s\n{y/E\ l v))LpZ(u). 

n=even n=odd 

(35) 

The separation of this solution to even and odd terms, though in principle unnecessary, is crucial for 
our prescription for the CIS. As stated before, this solution is valid only for small v. It is obvious that 
the presence of the odd terms of v dose not have any effect on the form of the initial wave function but 
they are responsible for the slope of the wave function at V = 0, and vice versa for the even terms. The 
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general initial conditions can now be written as 

tf(«,0) - ^A„^( w ) + ^C n ^H (36) 

even odd 

^(u,0) = iJ2 B nVE^Vl(u)+iJ2 Dn ^Vn(u), (37) 

eten odd 

where prime denotes the derivative with respect to v. Obviously a complete description of the problem 
would include the specification of both these quantities. However, given only the initial condition on the 
wave function, we show there is a CIS which produces a canonical wave packet with all the aforementioned 
desired properties. We can qualitatively describe the prescription for this case as setting the functional 
form of the odd undetermined coefficients to be the same as the even determined ones and vice versa. 
This means that the coefficients that determine CIS i.e. B n for n even and D n for n odd, are chosen as 
[22] 

B n = C n for n even D n = A n for n odd (38) 



In other word, by specifying the initial wave function, the prescription (]38[) automatically construct the 
appropriate initial slope which coincide well with the classical counterpart. Note that, although C n (A n ) 
is defined only for n odd (even), we can extend its definition to n even (odd) by choosing the same 
functional form. Now, using the canonical initial conditions (|36I37|) . we can determine As and construct 
the wave packet via equation (|29|) . 

The classical paths corresponding to these solutions can be obtained from (|16I17[) . The corresponding 
initial conditions for the classical case are 

u(0)=u , i>(0) = 0, u(0)=0, v(0)=v , (39) 

where the parameters uq and vq are adjusted so that the zero energy condition (|18|l is satisfied. 

For ease of comparison with FRW models [22] , we choose the same illustrative problem with io\ — 
ll>2 = u. Moreover, for all studied cases in this section (7 = 2, 4, 6) we choose the same coefficients as 

A n = e~3 |xl for n even B n = for n odd (40) 

V2 n nl 
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where \ is a free parameter. This choice of expansion coefficients obviously result in different initial 
conditions for various values of 7 (j3l))) . Note that, 7 = case is equivalent to FRW case with the constant 
spatial curvature [22] . 

For 7 = 2, the Lagrangian ([8|) can be written as 



which is equivalent to the flat FRW cosmological model, but with a modified cosmological constant, 
A' = A - 3/9. For this case the WDW equation JTHJ) reduces to 



This equation is in the form of an isotropic oscillator-ghost-oscillator and is separable in the configuration 
space variables. The general solution can thus be written as a sum over the product of simple harmonic 
oscillator wave functions with the same frequencies. The exact classical paths would be Lissajous figures 
in general. In particular, for io\ = u)2, by using the expansion coefficients in the form of equation (1401) . 
the corresponding classical paths are circles with radii \- The result is shown in the left part of Fig. [1] 
As can be seen from this figure, the parameters of the problem are chosen such that the initial state 
consists of two well separated peaks and this class of problems are the ones which are also amenable 
to a classical description. We should mention that there are a variety of different cases illustrated in 
Ref. [21j including u)\ =/= L02- Having precisely set the initial conditions for both the classical and quantum 
cosmology cases, we can now superimpose the results as illustrated in the right part of the Fig. [T] As 
can be seen from the figure, the classical-quantum correspondence is manifest. 

For 7 = 4, 6, the WDW equation (fT9| and the corresponding classical field equations (|16I17I18P are 
not exactly solvable. The classical equations can be solved numerically using customary algorithms, 
and the quantum cases using SM. In fact, for 7 > 2, the bound state solutions exist only for positive 
values of f} (j3 > 0). This is contrary to FRW case, where the bound state solutions can be obtained for 
positive, zero, and negative values of the spatial curvature |22) . Now, using the canonical prescription 
we can construct the wave packets which follow their counterpart classical trajectories. Figures 12131 show 




(41) 




(42) 
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Figure 1: 7 = 2 case: Left, the square of the wave packet \tp(u,v)\ 2 for w 2 = 1, j3 = 1, x = 3.5 and 
iV = 15. Right, the contour plot of the same figure with the classical path superimposed as the thick 
solid line. 

the resulting canonical wave packets and their classical trajectories for 7 = (4,6), and \ — (3-5,4), 
respectively. We have set (3=1, and used N = 15 basis functions to reproduce the wave packets. Note 
that, for these cases the parameter x corresponds to the classical initial position, but unlike the previous 
case, the classical pathes are no longer circles. 

An interesting feature of the Stephani Model is that it allows us to still have bound state solutions 
even with negative values of oj 2 . In fact, for 7 > 2 and (3 > 0, bound state solutions also exist for all 
positive values of A. Figure [4] shows the resulting classical and quantum mechanical solutions for 7 = 4, 
lo 2 = —1, and (3—1. Using the same expansion coefficients (|40[) . we have a slightly larger initial position 
description with respect to the previous case where 7 = 4, lo 2 — 1, and (3=1 (Fig. [2|). To be more 
specific, the classical initial positions for these cases are uq = 2.3, 2.4 for lo 2 = +1, —1, respectively. 

4 Causal interpretation 

To make the connection between the classical and quantum results more concrete, we can use de Broglie- 
Bohm interpretation of quantum mechanics. In this interpretation the wave function is written as 



*(w, v) = Re 



(43) 
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Figure 2:7 = 4 case: Left, the square of the wave packet \ip(u, v)\ 2 for uj 2 — 1, f3 = 1, x = 4 and iV = 15. 
Right, the contour plot of the same figure with the classical path superimposed as the thick solid line. 




Figure 3:7 = 6 case: Left, the square of the wave packet \ip(u, v)\ 2 for uj 2 = 1, (3 = 1, \ = 3 and N = 15. 
Right, the contour plot of the same figure with the classical path superimposed as the thick solid line. 
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Figure 4: 7 = 4 case: Left, the square of the wave packet \ifj(u,v)\ 2 for lu 2 = — 1, /? = 1, x = 4 and 
N = 15. Right, the contour plot of the same figure with the classical path superimposed as the thick 
solid line. 



where R = i?(w, v) and S = S(u, v) are real functions and satisfy the following equations 



i)S_ 

du 



1 2 2 



1 d 2 fl 



as 
dv 



2 2 1 



d 2 S 



^+I^ 2 - w2 ) (7+1)/3 

i? at;^ 4 



du 2 dv 2 du du dv dv 



0. 
0. 



(44) 
(45) 



To write R and S, it is more appropriate to separate the real and imaginary parts of the wave packet 



^{u, v) = x(u, v) + iy(u, v), 



(46) 



where x, y are real functions of u and v. Using (|4"5| we have 



R — \/ x 2 + y 2 



y 

S = arctan( — ). 

x 



On the other hand, the Bohmian trajectories are governed by 

dS_ 

du ' 
dS 



(47) 
(48) 



Pu — o ; 

du 



Pv 



dv 



(49) 
(50) 



where p u and p v are the momenta conjugate to u and v variables, respectively. Therefore, the Hamil- 
tonian constraint (H — 0) is again satisfied, but in the presence of the modified potential (|4"4"1) . The 
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Bohmian equations of motion take the form 

u = ^T7^i(!)' (51) 

where a; and y are known functions of u and v p5|). These differential equations can be solved numerically 
to find the time evolution of u and v. 

Using the explicit form of the wave packets, these differential equations can be solved numerically 
to find the time evolution of u and v. First, consider the case when 7 = 2. In this cases, it is apparent 
that the full potential for u (v) is no longer equal to u 2 (v 2 ) but is u 2 — j^t^t {v 2 — jzjfjr)- hi the 
right part of figures O and [51 we have shown the classical and Bohmian trajectories together for two 
different choices of initial wave function (A(n) — | e~ x2 / 4 , A(n) — J^* ; e~ x2 / 4 ). We see that 

the Bohmian trajectories are in good agreement with the classical counterparts. Now, let us find the 
quantum potential for instance in u direction along the Bohmian trajectories which is given by 

\8 2 R_ x' 2 + xx"+y' 2 + yy" f xx' + yy' \ 2 
Q R 8u 2 x 2 +y 2 \ x 2 +y 2 ) [ ' 

where prime denotes the derivative with respect to u. Figure [7] shows the classical (Vc) and quantum 
(Vq) potentials for two mentioned initial conditions. In particular, for A(n) — —^==e~ x / 4 , we found 
that for x ^ 3 (where \ is also the classical radius of motion for this choice of expansion coefficients) 
the functional form of the quantum potential is Vq — Vq(x/x) with the maximum value at x = %. This 
means that 

T^«-2> fOT X>3. (54) 

Moreover, as indicated in Fig. [SJ initial Bohmian velocity coincides well with the classical counterpart 
for large x which is compatible with the smallness of the quantum potential (|54[) . In fact, for this choice 
of expansion coefficients, the initial wave function consists of two lumps centered at u = ±x as (Fig- E]) 

°) = 1^74 (V (U - X)2/2 + e-^ 2 ' 2 ) . (55) 
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U 

Figure 5: 7 = 2 case: Left, the square of the wave packet \^p(u,v)\ 2 for A(n) = j e~ x / 4 and x = 5. 
Right, the classical (dashed line) and Bohmian (solid line) trajectories. 

Therefore, the complete classical and quantum correspondence occurs when there is no significant overlap 
between the two pieces of ift(u, 0). This means that to have a good correspondence for small radii, we 
need to choose a different set of coefficients or initial wave function which leads to a more localized wave 
function with infinitesimal overlap between its parts. We can also use causal interpretation for other 
cases. In particular, Fig. 1101 shows the obtained Bohmian positions versus time (i.e. u(t)) for 7 = 4 and 
7 = 6 which coincide well with their classical counterparts. 

5 Conclusions 

We have described a Stephani type cosmology near its symmetry center leading to classical dynamical 
equations given by (|16M18[) and the corresponding WDW equation represented by (TT!)|) . All these equa- 
tions are not exactly solvable and we have solved these equations numerically by an implementation of 
the SM for the quantum cosmology cases. We then constructed the wave packets via canonical proposal 
which exhibit a good classical-quantum correspondence. This method proposes a particular connection 
between position and momentum distributions which correspond to their classical quantities and respect 
to the uncertainty principle at the same time. Here, using canonical prescription, we tried to construct 
the wave packets which peak around the classical trajectories and simulate their classical counterparts. 
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Figure 6: 7 = 2 case: Left, the square of the wave packet \if)(u, v)\ 2 for A(n) — ™£ n] e x l A and x = 4. 
Right, the classical (clashed line) and Bohmian (solid line) trajectories. 
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Figure 7: fi(u,v): Classical (solid line) and quantum mechanical (dashed line) potentials for two types 
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Figure 8: 7 = 2 case: Classical (dashed line) and Bohmian (solid line) initial velocity versus x f° r 
A( n ) = -^L=e-* 2 / 4 . 
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Figure 9: 7 = 2 case: Left, the initial wave function ip(u,0) and right, the initial slope of the wave 



function iip'(u, 0) for A(n) 



= e x / 4 and x = 3.5. 
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We have also studied the situation using de-Broglie Bohm interpretation of quantum mechanics to quan- 
tify our purpose of classical and quantum correspondence and showed that the Bohmian positions and 
momenta coincide well with their classical values upon choosing arbitrary but appropriate initial condi- 
tions. Moreover, We showed that, in some cases, contrary to FRW cases, the bound state solutions also 
exist for all positive values of the cosmological constant. 
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